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Abstract. We apply a principal component analysis 
(PCA) to the spectra of each of the 18 Seyfert 1-like ob- 
jects observed more than 15 times by the international ul- 
traviolet explorer (lUE) from 1978 until the end of 1991. 
PCA allows us to decompose the Lya and Civ A1549 
emission-lines of each object into two components with 
uncorrelated variations. We find that the principal com- 
ponent describes correlated continuum and line variations, 
whereas the rest component reveals the line-parts that do 
not vary in tune with the continuum. A cross-correlation 
analysis in NGC 5548 reveals that variable line features 
in the rest component follow the continuum by about 25 
days. The symmetry usually observed in the component's 
line profiles excludes that the velocity field of the line- 
emitting region is dominated by radial motion. In some 
objects, the principal component reveals a clearly double- 
peaked line profile. 

The results suggest that the variable line feature in the 
principal and the rest component is respectively emitted 
in the inner and outer parts of an extended broad-line re- 
gion (BLR). The fraction of the line that varies together 
with the continuum seems to decrease significantly with 
increasing luminosity. We interpret this trend as a conse- 
quence of a larger BLR in more luminous objects. Finally, 
we propose a solution to the C iv/Lya ratio problem with- 
out referring to optically thin emission-line clouds. 

Key words: Line: profiles - Galaxies: active - quasars: 
emission lines - Galaxies: Seyfert - Ultraviolet: galaxies 



1. Introduction 

Intensive studies of some Seyfert 1 galaxies have shown 
that emission-line variations are clearly a response to con- 
tinuum variations (e.g. NGC 5548: Clavel et al. 1991; NGC 
3783: Reichert et al. 1994; Fairall 9: Rodrfguez-Pascual et 
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al. 1997). This gives strong support to the standard idea 
that emission-lines in active galactic nuclei (AGN) arise 
from photoionization of many small dense clouds moving 
around a central ionizing source (Netzer 1990). The nature 
of these clouds is still unclear, but an attractive possibil- 
ity is that they might be the envelopes of bloated stars 
(Alexander 1997 and references therein). 

To understand the structure and the kinematics of 
the emission-line region, it is of great interest to be able 
to decompose the emission-lines into several independent 
components. Most previous attempts decomposed the line 
profile into several Gaussian components (e.g. Wamsteker 
et al. 1990). Such decompositions do not ensure that the 
components are really varying independently. 

An interesting way to overcome this problem is to de- 
compose the emission-lines by a principal component anal- 
ysis (PCA). This was first done by Mittaz et al. (1990), 
who applied a PCA to the lUE spectra of NGC 4151. They 
concluded that PCA has only a fairly limited application 
to the analysis of spectral variability, because they could 
not find the physical meaning of most components. The 
analysis we propose here is similar except that we con- 
sider only two components: the principal component and 
the rest component, which comprises all minor variations. 
Doing so, the physical interpretation of the components 
becomes much easier and PCA turns out to be a promis- 
ing way to study emission-line variations and hence the 
structure and the kinematics of the line emitting region. 

The great advantage of the PCA, as compared to re- 
verberation mapping (Peterson 1993) or cross-correlation 
methods, is that it does not require very well-sampled ob- 
servations. PCA allows us to study objects that have much 
too sparse observations to perform a meaningful cross- 
correlation between continuum and line fluxes. Our sam- 
ple of 18 objects is to our knowledge the greatest on which 
line profile variations were studied in a consistent way. It 
gives us the opportunity to identify common behaviors 
among the AGN diversity and to look for possible trends 
with the object luminosity. 
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2. Observational material 

All spectra were down-loaded from the uniform low dis- 
persion archive (ULDA) version 4.0 in the Swiss ULDA 
national host in Lausanne. This database contains al- 
most all the low dispersion (resolution ~ 6 A) ultraviolet 
(UV) spectra obtained by the lUE satellite from 1978 to 
the end of 1991 (Courvoisier & Pahani 1992). The spec- 
tra of the ULDA database are processed with the lUE- 
SIPS reduction software described by Wamsteker et al. 
(1989) and references therein. We arc aware that the su- 
perior NEWSIPS reduction software would improve the 
signal-to-noise ratio of the spectra. However, the lUE fi- 
nal archive (lUEFA) docs not yet contain all the spectra 
used in this work and we considered that the improvement 
was not worth repeating the reduction for the more than 
thousand spectra used in this work. 

We considered the 18 Seyfert 1-like objects (quasars 
and Seyfert galaxies) with more than 15 good quality spec- 
tra observed through the large aperture (10" x 20") of the 
short wavelength prime (SWP) camera. Each spectrum 
was viewed separately and we excluded sky and noisy spec- 
tra, as well as some spectra with bright spots. The redshift 
of 3C 273 is so that the C iv A1549 line coincides with a 
reseau mark at Aobs 1792-1796 A, its most spoiled spectra 
were also excluded. 

Because of the lUE wavelength scale uncertainty due 
to small displacements of the target within the lUE aper- 
ture, we realigned all spectra according to the peak of 
the geocoronal Lya A1216 emission fitted by a Gaussian. 
This was done for all objects except NGC 4151, because 
the geocoronal emission is contaminated by its small red- 
shifted (z = 0.0033) Lya line. We choose the geocoronal 
emission rather than emission-lines from the object to 
avoid problems with double-peaked profiles and to use 
an identical procedure for all objects. The dispersion of 
the shifts that we measured is ~1.5A, but shifts up to 
~ 5 A are observed. This is enough to introduce spurious 
antisymmetric; line profile variations (Trier & Courvoisier 
1997). Finally, the spectra were rebinned into 1 A bins for 
computational convenience. Throughout this paper, the 
flux density is expressed in 10~"'^^ergs^^ cm^^ A and 
the integrated line flux is expressed in 10"^'' ergs~^ cm~^. 



3. Principal component analysis 

PCA is a mathematical tool, which can reduce a multi- 
dimensional data set into a small number of linearly in- 
dependent variables. A general description of the PCA 
method can be found in Kendall & Stuart (1976) or JoUifFe 
(1986). The special case of PCA applied to AGN spectra 
is described by Mittaz et al. (1990) or by Francis et al. 
(1992) for object-to-object variations. 
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Fig. 1. Principal component analysis of the CiV A1549 line 
in NGC 3783. a decomposition of the mean spectrum (points) 
into the principal component (solid line) and the rest com- 
ponent (dotted line), b similar decomposition of the variance 
spectrum. The flux is expressed in lO"'^* ergs~^ cm~^ A 



3.1. Description of the method 

The analysis we applied to the spectra of each object was 
done as follows. We constructed a m x n matrix con- 
taining the fluxes of the m observations in the n 1 A wave- 
length bins Ai {i = 1 . . .n) between Aobs 1229 and 1948 A. 
These 720 bins represent nearly the entire SWP range 
except the geocoronal Lya emission bellow 1229 A. The 
fluxes in each Aj form a base {/(Aj)} of a n-dimensional 
space, in which the m observations are represented by m 
points. 

The aim of the PCA is to form a new base {/(A^}, 
in which only a few number of the n vectors /(AQ are 
relevant to describe the variations. {/(A-)} is the base in 
which the expression of the covariance matrix C of is 
diagonal. The vectors /(AQ are thus the eigenvectors of C. 
They are called components, since they describe linearly 
independent forms of variations and can be represented as 
spectra (each component is a linear combination of the n 
initial /(A^)). For a given component k, the eigenvalue of 
C is equal to the flux variance (T^ in this component. The 
contribution of the component k to the total variability 
is given by its relative importance Ik defined by: Ik = 
(Tk/ J27=i (^i=^ fc/Trace(C). The principal component is 
thus the component with the highest Ik and its direction 
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in the n-space follows the maximal elongation of the m 

observation points. 

The matrix J-"', which is the expression of J- in the new 
base {/(A-)}, contains the fluxes of the m observations in 
the n "pseudo- wavelengths" A^. Wc therefore have the n 
lightcurves of the components from which we derive the 
mean flux ( / ) j and the flux variance (T^ in each compo- 
nent i. The raw components arc normalized to unity and 
their sign is undefined. In order to be able to add or sub- 
tract components, we have to normalize them properly. 
Two possible additive normalizations are (a) to multiply 
the components by their mean flux {f)i or (b) to multi- 
ply the squared components by their flux variance (7?. The 
total component (i.e. the sum of all components) is then 
equal to (a) the mean spectrum or (b) the variance spec- 
trum (i.e. respectively the mean flux and the flux variance 
in each Aj). 

In our analysis, we concentrate on only two compo- 
nents: the principal component and the rest component, 
which is constructed by subtracting the principal from the 
total component and is therefore the sum of the n — 1 mi- 
nor components. In Fig. 1, we show the decomposition into 
these two components of (a) the mean and (b) the vari- 
ance spectrum taken around the C IV A1549 line of NGC 
3783. The lightcurve of the total component is obtained 
by a vectorial addition of the lightcurves in the n com- 
ponents (i.e. for each observation the squared flux in the 
total component is the sum of the squared fluxes in the n 
components). 

3.2. Properties of the components 

The power of the PCA lies in the diagonalization of the co- 
variance matrix C of ^. A diagonal covariance matrix en- 
sures that all structures constituting a single PCA compo- 
nent vary simultaneously, while the different components 
have variations that are completely uncorrelated at zero 
lag. 

We applied the PCA on the widest possible spectral 
range to include a significant fraction of the continuum in 
the analysis. This ensures that the principal component re- 
flects the main continuum variations, since the continuum 
is known to vary with larger amplitudes than the emission- 
lines. Figure lb illustrates that in NGC 3783 nearly all 
continuum variations are described by the principal com- 
ponent alone. Under this assumption, we can physically 
interpret the line profiles in the two components. The line 
profile in the principal component shows which line-part 
varies simultaneously with the continuum, whereas the 
line profile in the rest component shows line-parts that 
do not vary in tune with the continuum. 

To test this interpretation and to further investigate 
the relationship between the continuum and our two com- 
ponents, we applied the PCA method to only the well 
sampled observations of the 1989 monitoring campaign 
on NGC 5548 (Clavel et al. 1991). In Fig. 2a, we compare 
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Fig. 2. PCA study of the 1989 monitoring campaign on 
NGC 5548. a logarithmic lightcurves of the continuum flux 
at A rest 1350 A from Clavel et al. (1991) (star points), of the 
principal component (filled points) and of the rest component 
(open points) . b cross-correlations of the principal component 
(filled points) and the rest component (open points) with the 
continuum lightcurve at A rest 1350 A. c PCA decomposition of 
the C IV line as in Fig. la 



the principal and the rest component lightcurves with the 
original continuum lightcurve at A rest 1350 A from Clavel 

et al. (1991) in units of 10~^^ erg s~^ cm~^ A . 

We analysed the correlation of the two components 
with the UV continuum using the interpolated cross- 
correlation function (ICCF) of Gaskell & Sparke (1986) 
as described by White & Peterson (1994). The results 
show that the principal component variations are indeed 
strongly correlated with those of the continuum flux at 
zero lag, while the rest component variations are com- 
pletely uncorrelated at zero lag with those of the con- 
tinuum. However, the cross-correlation between the rest 
component and the UV continuum suggests that the vari- 
able line-parts in the rest component follow the continuum 
with a delay of the order of 25 days. This lag is twice as 
large as the lag between the continuum and the integrated 
Lya and Civ lines (8-16 days) (Clavel et al. 1991), which 
means that the PCA divides the line into two parts: a line- 
part that responds to the continuum with a very small de- 
lay { ^5 days) and a line-part that responds with a much 
longer delay (~25 days). 
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The Civ line profile of these two parts is shown 
in Fig. 2c. The full width at half maximum (FWHM) 
of the line-part in the rest component is much nar- 
rower (1850kms~^) than in the principal component 
(7300kms-i). However, the full width at 1/4 maximum 
of the line-part in the rest component (6 650 km s~ ^ ) shows 
that high velocities are also present. 

This analysis suggests that the line-part in the prin- 
cipal component is emitted in the inner broad-line region 
(BLR), since it responds with large amplitudes and nearly 
zero lag to the continuum, whereas the line-part in the 
rest component is both constituted of a less varying con- 
tribution from outer parts of the BLR that respond with 
a greater delay to the continuum and of a non-varying 
contribution from the narrow- line region (NLR). 

The lag of ~ 25 days found for the rest component in 
NGC 5548 is comparable to the lag of the Cm] A1909 
line (26-32 days) (Clavel et al. 1991) and to the Balmer- 
line lags (^20 days for H/3) (e.g. Peterson et al. 1991). 
It is therefore possible that the C iv line-part in the rest 
component is emitted at about the same distance from the 
ionizing continuum than the C III ] or the H/3 line. 

In other objects, the temporal sampling of the obser- 
vations does not allow us to determine a clear lag between 
the two components and the continuum. However, this will 
be possible with the other monitoring campaign data, that 
will soon be available in the lUEFA. 

4. Definition of the variability parameters 

From the PCA results, it is possible to derive line and 
continuum variability parameters. The mean continuum 
if \,cont)k and line {fiine)k fluxes are derived from the 
component k normalized by its mean flux ( / )fc (Fig. la). 
The continuum (J A.cont and line (7 unc flux dispersions arc 
estimated from the principal component normalized by its 
flux dispersion (Jprin- This derivation is motivated by the 
theoretical argument exposed in Sect. 4.1. The determina- 
tion of the continuum and the method of line integration 
is described in Sect. 4.2. 

4-1- Theoretical basis 

In the following, the temporal average (./(A)) and the 
temporal dispersion (T(A) of the n spectra f{X,ti) mea- 
sured at time ti (i = 1 . . .n) are defined as 

1 " 

(/(A)) = -Y.f{\U) and (1) 

1=1 

1 " 

^'(A) = -—r^{Af{X,U))\ (2) 
" ^ i=i 

where A/(A, U) = /(A, U) - { /(A) ). At each U, the total 
spectrum is the sum of a continuum and a line contribu- 
tion /tot(A,tj) = /coiit(A,ii) + /iine(A,ti). This relation 



combined with Eqs. (1) & (2) gives 

(/tot (A)) = (/cont(A)) + (/line (A)) (3) 
O-tot(A) = CrlontW + line (A) 

2 " 

+ ——r A/cont(A, U) ■ A/line(A, U) (4) 

Eq. (4) shows that contrary to the mean spectrum, the 
variance spectrum can not be simply decomposed into a 
continuum and a line contribution. 

It would however be possible to decompose the disper- 
sion spectrum (i.e. the flux dispersion in each Aj) into a 
continuum and a line contribution if 

A/line(A, U) = R{X) ■ A/eont(A, U) , (5) 

because this would imply that (7 linc(A) = -R(A) • Crcont(A) 
via Eq. (2) and thus that Eq. (4) reduces to 

^^tot(A) = ^J(Jl^^dX) + O-Le(A) + 2R{X)al,^^{X) 

= crcont(A) + crii„e(A). (6) 

Equation (5) would correspond to substitute a time 
independent response function R{X) for actual convolution 
with the transfer function. Since there is a signiflcant lag 
between the line and the continuum variations, Eq. (5) is 
incorrect and thus we can not derive the line variability 
(7 line from the dispersion spectrum via Eq. (6). 

However, we have shown in Sect. 3.2 that the princi- 
pal component is correlated with the continuum at zero 
lag. This is not only verifled in NGC 5548, but is a direct 
consequence of the way we applied the PCA. Therefore, 
Eq. (5) is correct for the line-part in the principal com- 
ponent and thus, according to Eq. (6), we can estimate 
the line variability <7 une by integrating the line-part in 
the principal component normalized by its flux dispersion 
Cprin- This estimation gives a lower limit to the actual 
line variability, since it does not consider the delayed line 
variations, but has the advantage of excluding spurious 
variability due to noisy spectra. 

4-2. Continuum and emission-line determination 

On the basis of Eqs. (3) & (6), we can derive both the 
mean and the dispersion of the line flux by integrating the 
emission-line above a deflned continuum in the appropri- 
ate component. To do this, we first defined five 30 A bands 
that are usually free of emission or absorption features in 
the rest frame of the source at Arest 1120-1150 A, 1320- 
1350 A, 1430-1460A, 1700-1730 A, 1810-1840 A. The two 
last bands are both outside of the component domain 
Aobs 1229-1948 A in 3C 273 and GQ Comae, due to their 
higher redshift. For these two objects, we had to define a 
new band from Arest 1665 A to the end of the component, 
despite the possible He ii A1640 and O iii ] A1663 contam- 
ination at those wavelengths. Doing so, there are always 
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Table 1. List of the objects in our sample with some characteristic parameters derived from the PCA. log(L) is the logarithm 
of the luminosity expressed in ergs~^ A ^, the FWHM is expressed in kms~^. Other quantities have no units. The parameters 
are defined in Sect. 5.1 
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GQ Comae 


17 


0.1653 


41.93 


0.66 


0.18 


0.40±0.05 


0.38±0.11 


4 750 


5 000 


3 100 


2 650 


Mrk 279 


23 


0.0297 
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0.11 
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3 700 


7 000 
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37 


0.0355 


41.60 
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0.16 
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0.53±0.03 


6650 


5 350 


3 300 


2 400 


Mrk 926 


18 


0.048 


41.54 


0.85 


0.14 


0.43±0.03 


0.61±0.02 


7600 


10 650 


5 000 


5 050 


NCC 3516 


30 


0.009 


39.73 


0.81 


0.18 


0.71±0.02 


0.60±0.01 


6100 


8 900 


2 700 


2 300 


NGC 3783 


55 


0.0096 


40.24 


0.87 


0.18 


0.97±0.03 


0.60±0.01 


7 400 


4 650 


2 300 


2 200 


NCC 4151 


314 


0.0033 


39.60 


0.89 


0.13 


0.94±0.08 


0.82±0.01 


8 000 


7 650 




1850 


NGC 4593 


25 


0.0087 


39.58 


0.76 


0.17 


0.79±0.04 


0.68±0.03 


10 300 


9 250 


3 800 


3 500 


NGC 5548 


154 


0.017 


40.71 


0.78 


0.17 


0.91±0.01 


0.63±0.01 


6 900 


7 450 


1500 


1850 


NGC 7469 


20 


0.0165 


40.77 


0.51 


0.42 


0.27±0.04 


0.33±0.05 


7 000 


7 500 


2 000 


2 850 



four bands that constrain the continuuin, whatever the 
redshift of the object is. 

The continuum is estimated by fitting a straight line 
through the points in these four bands. Wc used the or- 
dinary least-squares regression "OLS(X|Y)" of Isobe et 
al. (1990) with its uncertainties on the slope and on the 
intercept. To take into account systematic errors, the un- 
certainty on the continuum is assumed to be twice the 
uncertainty on the fit that was calculated, as all uncer- 
tainties in this paper, according to the general formula 




where / is a function of the n variables Xi with uncertain- 
ties (7 a;.. 

The Lya and C IV emission-lines are integrated above 
the continuum in the velocity ranges [ — 15 000 ; -|- 15 000 ] 
and [ —15 000 ; + 12 000 ] respectively, as far as enabled by 
the redshift of the object. We choose these ranges in order 
to include the Nv A 1240 line in the Lya A1216 integra- 
tion and to exclude the Heii A1640 and O ill] A1663 lines 
from the C IV A1549 integration. In the following, the Lya 
line actually refers to the Lya -|- N V line. The uncertainty 
on the line contribution is determined by integrating the 
line above the fitted continuum plus and minus its uncer- 
tainty. The continuum parameters are arbitrary defined at 
A rest 1350 A. 



5. Results 

In the appendix, wc show the PCA decomposition of the 
Civ and Lya lines for all objects. If the component's line 
profiles are similar for the C iv and the Lya line of a same 
object, they can be very different from one AGN to the 
other. This great diversity forced us to consider only very 
simple parameters to quantify the line properties of the 
components in order to find significant trends. The general 
results and possible trends with luminosity are presented 
in Sect. 5.1, whereas comments on individual objects are 
given in Sect. 5.2. 

5.1. Global trends 

Tables 1 & 2 display the most meaningful quantities de- 
rived from our analysis together with the number n of con- 
sidered spectra and the redshift z of the object as given 
by the set of identifications, measurements and bibliogra- 
phy for astronomical data (SIMBAD). log(L) is the loga- 
rithm of the luminosity density at A rest 1350 A expressed 
m ergs" -iA-\ It was derived from the fitted continuum 
assuming a Hubble constant of Hq = 50kms~^ Mpc~^. 

To quantify the variability in the principal component 
with respect to other variations, we choose two param- 
eters: /prin and (7 rest /C prin, which are respectively the 
relative importance of the principal component and the 
ratio of the fiux dispersions in the rest and in the princi- 
pal component, as defined in Sect. 3.1. Both parameters 
show that in general the dominant variations are well de- 
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Fig. 3. Fraction of the emission-line present in the principal 
component as a function of the object luminosity. Open and 
filled points represent respectively Lya and C iv linos. The two 
dotted curves represent the 3-a uncertainty on the bisector 
linear regression of Isobe et al. (1990) 



Fig. 4. Distribution of the FWHM of the emission-lines in the 
rest component (a) and the principal component (b). The CiV 
lines arc shaded or twice shaded if they are doublo-peakcd. 
The mean (dashed line) and the dispersion (dotted lines) of 
the distribution are shown 



scribed by the principal component alone (/prin ~ 50 % 
and (J rest /cr prin ^ 50%). 

According to the discussion in Sect. 3.2, the fraction 
Tj of the line that varies together with the continuum is 
given by the ratio of the mean line fluxes in the principal 
and in the total component f] = { fune ) prin/ ( /line ) tot - Its 
dependence with luminosity shown in Fig. 3 has a slope of 
— 0.28 ± 0.06/decade. Here, as in the following, the slope 
was determined by using the "OLS bisector" unweighted 
linear regression of Isobe et al. (1990), as advised by these 
authors and the trend is highlighted in the figure by the 
3-a uncertainty on the regression. We applied a Spear- 
man's test (e.g. Bevington 1969) to determine whether the 
correlation suggested by Fig. 3 is significant or not. The 
probability that such a correlation could occur by chance 
is 1.6%. 

We quantified the width of the Lya and the C iv lines 

by evaluating their FWHM both in the principal and in 
the rest component. The FWHM values displayed in Ta- 
ble 1 are expressed in kms~^ with an uncertainty on the 
measure of ±100 km s"-*^. Some noisy components had to 
be smoothed to determine properly their width and in low 
redshift objects, for which the blue side of the Lya line is 
not in the component, the FWHM was extrapolated from 



the half width at half maximum. Even so, we could not 
determine a meaningful width for the Lya line in the rest 
component of NGC 4151, because of geocoronal Lya con- 
tamination and for the C IV line in the principal compo- 
nent of 3C 273, because there is nearly no line feature in 
this component. 

The distribution of the FWHM in Fig. 4 shows that 
there is no obvious difference between the CiV and the 
Lya line. The average width of the two lines is twice as 
large in the principal component (6 850 km s~^) as in the 
rest component (3 400kms~^). Individually, the line is al- 
ways broader in the principal component than in the rest 
component, except for the Lya line in 3C 120 and the C IV 
line in Mrk 279. This clear result shows that the line-part 
that varies with the continuum is broader in general than 
the line-part that varies less. In many objects however, 
the line in the rest component is broader than the typi- 
cal width (300-1 OOOkms-i) of a line emitted in the NLR 
(Netzer 1990). 

The Civ line profile in the principal component is 

clearly double-peaked in 3C 390.3, NGC 3516 and NGC 
4151. A double-peaked profile is known to be the signature 
of a thin rotating disk viewed close to edge-on (Welsh & 
Home 1991). Recently, Goad & Wanders (1996) showed 
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Table 2. The relative variability of the continuum measured 



at A rest 1350 A and of the integrated Lya and C IV lines 
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that double-peaked profiles can also originate due to a 
non-uniform lighting of a spherical BLR by a predomi- 
nantly biconical continuum emission. Their model is able 
to reproduce a wide range of observed profiles and predicts 
that the FWHM should be generally larger in double- than 
in single-peaked profiles. The FWHM that wo measured 
for our three CiV double-peaked profiles (^^ 8 000 km s^^) 
are among the highest in the sample, in good agreement 
with their prediction. The line profile in Mrk 509 can be 
seen as the transition between double- and single-peaked 
profiles. 

We quantified the variability of the continuum and of 

the Lya and C IV lines by the ratio of the flux disper- 
sion on the mean flux (7/{f). This ratio corresponds to 
the parameter Fvar first defined by Clavel et al. (1991) 
to describe the flux variability. We calculated the values 
shown in Table 2 as explained in Sect. 4. The luminosity 
dependence of the line variability shown in Fig. 5a has a 
slope of —0.16 dz 0.03/dccadc. The Spearman's test prob- 
ability that such an anticorrelation could occur by chance 
is 1.8%, but rises to 10%, when we consider only the 
Lya or the C IV line and reaches 70 % for the continuum 
variability. This absence of correlation between the con- 
tinuum variability and the luminosity is in contradiction 
with Paltani & Courvoisier (1994), who found a strongly 
significant anticorrelation in a sample of 72 Seyfert-like 
objects with a slope of — 0.046/decade. The fact that this 
trend remains hidden to us with a sub-sample of 18 ob- 
jects is most probably due to its weak slope and to the high 
scatter of the points. This increases our confidence in the 
significant trends that we found for the two emission-lines. 



Fig. 5. Luminosity dependence of the line variability (a) and a 
parameter describing the line-to-continuum response (b). Open 
and filled points represent respectively Lya and IV lines. The 
dotted curves are defined as in Fig. 3 



The ratio of the line variability on the continuum vari- 
ability (Cr/(/))ii„e/(Cr/(/))i35o is perhaps even more 
meaningful, since it describes the linc-to-continuum re- 
sponse. Its luminosity dependence shown in Fig. 5b has a 
steeper slope (— 0.24± 0.05/dccadc) than that of the line 
variability, but the Spearman's test probability is similar 
(1.7%) and the scatter of the points is higher. The great 
similarity with the relation in Fig. 3 is due to the fact 
that T] and {(J/{f)) unc/ {(7 /{,f ))i35o are closely related, 
as shown by the following theoretical argument. 

Since the PCA results suggest that a fraction Tj of the 
line varies in tune with the continuum, while the remain- 
ing fraction (1 — Tj) is nearly constant, we can write 

/li„e(fi) ~ i?[r//A,co„t(t.) + (l-^)(/A,co„t)] (8) 

-cont )+77A/A,co„t(ii)], (9) 

where R\s a, response constant and Af{ti) = f{ti) — ( / ). 
We can now equal the first and the second term in Eq. (9) 
with (/line) and AfuneiU) respectively: 



( /line ) 
A/li„e(ti) 



R ( /a, cont ) 
Rr]Afx,cont{ti)). 



(10) 

(11) 
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In accordance with Eq. (2), Eq. (11) implies that (7 une 
RT](J A,cont, which together with Eq. (10) leads to 



7] 



{(7/{f)) 



line 



((J/(/))A,cont ■ 



(12) 



This straightforward calculation shows that both Tj and 
('-''/(/)) iine/(C/(/))i35o describe in a way the line-to- 
continuum response. 

In all but two objects (3C 390.3 and Mrk 926), the Hne 
response is stronger within the uncertainties for the Lya 
than for the C IV line. This naturally leads to a decrease of 
the Civ/Lya ratio with increasing continuum flux, as it 
was observed in some Seyfert galaxies (Peterson 1993). 3C 
390.3 was until now the only object in which the C iv/Lya 
ratio was observed to increase with increasing continuum 
flux (Wamsteker et al. 1997). Our results confirm this pe- 
culiar behavior in 3C 390.3 and predict a similar one in 
Mrk 926. 

We did not quantify line profile asymmetries, but a 
look at the components displayed in the appendix shows 
that most of their line profiles arc roughly symmetric. The 
observed absence of strong line asymmetries in the com- 
ponents of most objects cxchidcs that their line-emitting 
region is dominated by iiifall or outflow (Rosenblatt et 
al. 1994). However, minor line asymmetries might reflect 
some radial motion within a predominantly chaotic or ro- 
tational velocity field. 

5.2. Comments on individual objects 

Apart from 3C 390.3, which has a different C iv/Lya ratio 

behavior than most AGN (Sect. 5.1), other objects have 
some peculiarities. The most luminous radio- loud quasar 
in our sample is 3C 273, its principal component has only 
a very broad and nearly undetectable line feature. This 
leads to a very weak line variability, as was already no- 
ticed by Ulrich et al. (1993) for the Lya line. Since 3C 
273 is the only well enough observed luminous object, it 
is difficult to know whether or not its blazar characteristics 
are responsible for this peculiarity. However, the fact that 
3C 273 is well inside the obtained luminosity trend sug- 
gests that other luminous quasars would also have nearly 
constant emission-lines. 

Fairall 9, on the other hand, is an object which does 
not follow the general luminosity trend. Its line variabil- 
ity behavior is typical for a Seyfert galaxy, whereas its 
luminosity is that of a quasar. This was already pointed 
out by Rodrfguez-Pascual et al. (1997), who found that 
the emission-line lags in Fairall 9 (Lya: 14-20 days; Hell 
A1640: ;S 4 days) are comparable to those in NGC 5548, 
despite the difference in luminosity of a factor of ten. It 
suggests that the luminosity of Fairall 9 during the first 
years of lUE observations was extraordinary high and 
hence that its average luminosity is overestimated. 

3C 120 and NGC 7469 are two other objects that are 
always outside of the 3-(t uncertainty curves. Their line 



variability is very small, because they are the two only 
objects with a clear asymmetric line profile in their first 
component. The blue wing of the line is not correlated 
with the continuum in 3C 120, whereas it seems even an- 
ticorrelated in NGC 7469. Such asymmetries are quali- 
tatively consistent with an infalling BLR, since the blue 
wing is then emitted behind the continuum source having 
thus a higher emission-line lag than the red wing emitted 
in front of the source. However, both objects are not well 
enough observed to draw strong conclusions only based on 
these PCA results. 

6. Discussion 

Perhaps the most important result of the PCA is that 
a great fraction of the emission-line does not vary to- 
gether with the continuum. This fraction has a line profile 
that is usually too broad to originate only in the NLR 
and thus the classical division of the emission-line into 
a variable BLR and a constant NLR contribution has to 
be reviewed. The PCA decomposition gives support to 
the model proposed by Brotherton et al. (1994) in which 
the traditional BLR is divided in two components: a very 
broad line region (VBLR) and an intermediate line region 
(ILR), which have typical velocities of '~7000kms~^ and 
^ 2 OOOkms""'^, respectively. The line profile in the princi- 
pal component would then be the signature of the VBLR, 
whereas the rest component's profile would be the sig- 
nature of the ILR plus a constant contribution from the 
NLR. 

Under this assumption, the delay of ~25 days that 
we measured between the rest component and the UV 
continuum in NGC 5548 (Sect. 3.2) can be understood as 
the distance of the ILR. This result is in good agreement 
with Brotherton et al. (1994), who estimate that the ILR is 
about 10 times more distant from the ionizing continuum 
than the VBLR. We note however that the line profile in 
the rest component has wings which extend to significant 
fractions of those in the principal component. This implies 
the presence of high velocity gas in the ILR and hence 
suggests that the VBLR and the ILR are not completely 
separated. 

The ratio of the line on the continuum variability nat- 
urally decreases when the light-travel time across the line- 
emitting region t it is of the order of the typical variability 
timescale of the ionizing continuum tvar- If iit ^ ivar, 
the entire line-emitting region will respond as a whole to 
continuum variations, whereas if tit ^ ^var; different sub- 
regions will respond with different delays to the ionizing 
continuum averaging out the variations (Rosenblatt et al. 
1994). In our interpretation, the lines emitted in the ILR 
have weak variations, because the extended ILR implies 
that i It > t var , whereas the lines emitted in the VBLR vary 
in tune with the continuum, because this region is small 
enough to satisfy tit <tw A simple way to formally link 
the fraction Tj of the line in the principal component (i.e. 
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emitted in the VBLR) to the two characteristic timescales 
is to write 

7^~exp(-tit/ivar), (13) 

where t\i is now the Ught-crossing time of the BLR. Ae- 
cording to Eq. (13), we obtain with our data the foUowing 
relation with luminosity 

<itAvar=^-ln(77)ocLO-52±0-09^ (14) 

which shows that the BLR size increases with respect to 
ctvar as the luminosity increases. If we assume that t\t 
is proportional to the observed delays between the line 
and the continuum and thus that it varies roughly with 
the luminosity as tit oc L"-^ (Kaspi et al. 1996), the rela- 
tion in Eq. (14) suggests that the variability timescale of 
the ionizing continuum fyar is independent of the object 
luminosity. 

The usually observed decrease of the Civ/Lya ratio 
with increasing continuum flux was interpreted as being 
due to a population of optically thin clouds (Shields et al. 
1995). Such a model predicts that the broad emission-line 
component emitted in the inner BLR by optically thin 
clouds should be less variable than the narrower compo- 
nent emitted in the outer BLR by optically thick clouds. 
This prediction seems to be verified for the H/3 line in 
some Seyfcrt galaxies (Mrk 590: Peterson et al. 1993; Mrk 
335: Kassebaum et al. 1997), but was never observed to 
our knowledge for the Lya or the C IV line. Our result that 
the most varying line-part is usually also the broadest is 
in contradiction with such a model. 

We propose that the C iv/Lya ratio problem is due to 
a significant decrease of the variability timescale t^ai be- 
tween the ionizing thresholds of H I and C IV, respectively 
at 912 A (13.6 eV) and at 192 A (64.5 eV). If we assume 
that the light-crossing time of the BLR tn is similar for 
the two lines, this would lead to a better line-to-continuum 
response for the Lya than for the CiV line (Eq. (12) & 
(13)), implying a decrease of the Civ/Lya ratio with in- 
creasing continuum flux. 

The lack of far UV observations does not allow us to 
test this interpretation. However, there are evidences that 
the continuum varies more rapidly at higher energies than 
in the UV. Nandra et al. (1991) found that the X-ray flux 
in NGC 5548 varies significantly on timescales of hours 
and that it can vary by a factor of two in less than a few 
days. More recently, Marshall et al. (1997) found that the 
extreme UV (1 50-200 eV) can vary by a factor of two on 
timescales of 0.5 day. Such a decrease of the variability 
timescale toward shorter wavelength from the UV to the 
extreme UV domain is qualitatively in agreement with our 
interpretation. 
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Appendix A: PCA decomposition of the C iv line 



PCA decomposition of the mean C iv A1549 line profile (points) into the principal component (solid lino) and the rest component 
(dotted line). The abbreviated name of the object is shown in each panel. The zero velocity (dashed line) is defined by the 
redshift. On the flux axis, the first big tick is always at zero flux and small ticks are at every 10~^^ ergs~^ cm~^ A ^ . 
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Appendix B: PCA decomposition of the Lya line 

Same as in appendix A, but for the Lya line. The blue wing of the line could not always be analysed because of geocoronal Lya 
contamination. 
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